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Abstract 

We provide direct elementary proofs of several explicit expressions for Bernoulli 
numbers and Bernoulli polynomials. As a byproduct of our method of proof, 
we provide natural definitions for generalized Bernoulli numbers and poly- 
nomials of complex order. 
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1. Introduction 

As Gould pointed out in his survey article (HI, a well-known explicit for- 



mula for Bernoulli numbers, which dates back to Worpitzky |10J , is given by 
the double sum 



n— n -i— n \ J / 



(1) 

n=0 ' j=0 ^ ' ' 

where ( n ) is the binomial coefficient. The first few values of formula ([TJ are 
Bo = 15i = -1/2, B 2 = 1/6, B 3 = 0, B 4 = -1/30, B 5 = 0, B 6 = 1/42, 
B 7 = 0, B 8 = -1/30 etc. 

For the sake of convenience and to agree with our notation, the lower 
limits of summation in both sums in (TjQ) will be changed so that the above 
sum is given by the following equivalent form 
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n=l 7=1 VJ/ 



> 1. (2) 



.7 = 

If we define the forward differences A n (k) by 



An(*)=E(-l) i (")j*, (3) 

i=i 

equation (T5]) can be rewritten as 

£? fc = V— ^-A„(fc),£;>l. (4) 
z — f n + 1 

n=l 

In this note, our main result is proving the following two explicit formulas 
for Bernoulli numbers 

k 

B k = (-l) k+1 £,^ ^yA^),^!, (5) 

n=l ^ ' 

and 

fc+i , 

B k = ^ —A n (k + l),k> 0. (6) 



n 

n=l 



While formula ([I]) is a well-known explicit formula and several proofs 
have been provided by different authors, the two formulas (jSJ) and ([6]) are 
less-knowro 

As a byproduct of our method of proof we provide extensions to Bernoulli 
polynomials. We further provide natural definitions for generalized Bernoulli 
numbers and polynomials of complex order. 



1 Formula ([5]) has been given in 0, formula (37)]. The same formula is mentioned in 8, 
formula LXV on page 83] . Formula (|6]) is also mentioned in 0, formula LXIII on page 82] . 
The proofs in [§] use the identity A„(fc) = n(A n (k — 1) — A n _i(fc — 1)) and the property 
that odd-indexed Bernoulli numbers have zero values. In 0, page 83] the identity is written 
in terms of a\ as a k n = n(a^ _1 + a^z\) so that A n (fe) = (— l) n o* = (—l) n n\S(k, n), where 
S(k, n) are Stirling numbers of the second kind. 
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2. Fractional Derivatives 



Suppose that the function ip(t) is holomorphic and that lim t _ > ._ 00 if)(t) = 0. 
According to Laurent Q , the fractional derivative of order a e C between the 
points — oo and x E K of the function ^>(i) is given by the contour integral! 

where C is the Hankel contour consisting of the three parts C = C_ UC e U 
C + : a path which extends from (—00, — e), around the point x counter clock- 
wise on a circle of center the point x and of radius e and back to (— e, —00), 
where e is an arbitrarily small positive number. 

When Re(l + a) > 0, there is no ambiguity in the definition of I a i/j(x). 
The integrals along C_ and C + cancel each other. I a ip(x) is thus equal to 
the integral along C e and this integral can be easily evaluated by residue 
calculus. In particular, when a = 0, formula (J7|) is simply Cauchy's formula: 

J ^( x ) = ip( x ) : (8) 

and when a = n is a positive integer, Laurent's contour integral I a ijj(x) 
gives the classical derivative of ip(t) at the point x: 

I n ^( x )=^ n \ x ). (9) 

When Re(l + a) < 1, the portion of the integral along the circle C t is zero. 
The integral along the remaining portions of the contour is estimated using a 
proper determination of the multi- valued function (t — If we choose 

the cut long the semi-axis (— 00, x), then (t — x) -0 ^ 1 = e ~{ a + 1 )( lo &i. t ~ x )- l ' K ) 
along C- and (t - x)-*- 1 = e -("+i)(i°g(^)+^) a i ong C+, where log(t - x) is 
purely real when t — x > 0. Moreover, t = re~ t7T along C_ and t = re t7r along 
C + , as r varies from e to +00. The integral (J7J) becomes 

TV 1 -U n/\ /*°° 

I a ^(x) = 1 ; ( e ~ Q7ri - e°^) / ^(-r)(r - a;)^ 1 dr. (10) 



2m 

Finally, using the reflection formula of the Gamma function, we obtain 



2 Thc derivative can of course be denned for all x £ 
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I a i>(x) = — -/ ij){-r){r - x)~ a ~ l dr. (11) 

v ) 

3. Specializing to C( s ) 

For the particular = and a = — s, the fractional derivative of 

order — s at zero is given by 

/->(0) = ^-^ [ mt s ~ l dt (12) 
when Re(l — s) > 0, and by 



c 



I poo 

rSm = W)Jo ^~ t)f ' ldt (13) 

when Re(s) > 0. 

Now by an appropriate choice of ip(t), we will be able to write the Rie- 
mann zeta function as the fractional derivative of ip(t) at 0. Indeed, let 

^- t ) = ^ = s(^) = 7T^¥-T^- < 14 > 



We have shown in |3| that the Riemann zeta function has the integral 
representational 

1 f°° 

(s - l)C(s) = — - / ^(t)**- 1 dt, Re(s) > 0, (15) 

1 V s ) Jo 

and that for all s G C 

(s - l)C(a) = / mt 8 ' 1 dt, (16) 



where C is the same Hankel contour used in equation ([7]). Comparing 
with equations ( TT2l and (JE]), we easily obtain 

^-*W)U=(*-i)CW. (17) 



3 The definition of cf>(t) used here is different from the one in the cited paper. 
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That is, (s — 1)C( S ) is simply the fractional derivative of order — s of the 
function if)(t) at the origin. Furthermore, for an integer k > 2, the derivative 
of order k — 1 of ip(t) at the point t — is, by Laurent's definition, given by 

I'-'imu = ^ 1} (o) = v -^]mt~ k dt. (is) 

Having achieved this, we know also that the Bernoulli numbers are usually 
defined using the generating function 

l = £fr**' |t|<27r - (19) 

k=0 

These numbers can also be defined in terms of the function ip(t) instead, 
since we have 

k=l v ' 

Therefore, 



e*- 



B k = (-l)V (fc - 1) (0) = (-l) fc / fc - 1 [^W]* = o = (- 1 ) fc+lfc C(l-A:), or 
= _0(^)(o) = -/ fe - 1 [0(-t)] 4=o = A;C(l-A;). (21) 

This last equation is the basis of all our subsequent derivations. It relates 
the Bernoulli numbers, the Riemann-zeta function and fractional derivatives 
in a single equation. 



4. The First Explicit Formula for B k 



In [3|, we have also obtained a globally convergent series representation 
of the Riemann zeta function. It is given by the formula 

{s ~ l)C(s) = £ with (22) 

n=l 

n-l / _ i \ 

S n (s) = J2(-±) k [ , (A: + l)- s forn>2, (23) 

k=0 ^ ' 
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and Si(s) = 1. We have also shown that when Re(s) > the sum S n (s) 
can be rewritten as 

S n (s) = —--J {l-e-^e-Hr 1 ^, (24) 
and that the function </>(£) verifies 

n=l 

uniformljfl for < £ < oo . 

Because of equation ( 12"1"|) . the definition of fractional derivative (IT5|) and 
uniform convergence, we may interchange summation and integration inside 
the integral sig Thus, we may apply the fractional derivative operator 
termwise to obtain 



JEL T- S \(A — pt\n-Xpt] 

i- s [m] t=0 = i- s [<i>(-t)] t=0 = y: [[ w \ Jf =° - (27) 



71=1 

Particularly, for — s = k — 1 we obtain 



n=l 

71=1 

OO ^ ^ 

^n(n + l)^ [(1_e<)n]i = ' 



(28) 



71=1 



But ^[(1 - e*) n ] t=0 = if n > k + 1. Therefore, the infinite sum in (EE 
reduces to a finite sum 



4 Whcn Re(s) < 1, S n (s) can obviously be written as 



c 



2ni 

3 This have been rigourously proved in 



S n ( s ) = T (\ , s) / (1 - e^-Vt 8 - 1 dt. (26) 
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*" 1 d 



k 1 n (n\ d k 

^ n(n + l) (J ^ [e ' i]i =° 

n=l v ; j=0 VJ/ 

fc 1 " / \ 

n=l v y 7=0 yJ 7 



71=1 ' 

The last equation combined with ( |2T]) gives the explicit formula for 5*.: 

k 

B k = {-ir^ : £—K{k),k>l. (30) 

n=l ^ ' 

5. Bernoulli Polynomials -B fc (l — x) 

The generalization of the series formula ( 1221) for the Hurwitz zeta function 
£(s, x) defined for < x < 1 by 

C(s,x) = V- (31) 

v ' ^ (n - 1 + x) s v ' 

n=l v ' 

is given in j3| by the formula 



(s - l)C(s, x) = J2 Sn(s, x) (-— + , (32) 

n=i 

where ^(s,^) is the generalization of the sums S n (s): 



5 n (s,x) = ^(-1H ™ )(A; + x)- s forn > 2. (33) 

i — n \ 



fc=0 

There is also a corresponding integral given by 



(s - l)C(s, x) = — Mt)t s ~ l dt, (34) 

1 l s J </0 



where <j> x (t) is defined by 



te -xt e -{x-l)t ^ _ !) te -(x-l)t 

- ( et _ 1)2 - e * _ i + ^F31 
1 / -t \ £ - {x - 1)t + (x - ljte-C'-^ 



dt\e*-l/ e*-l 

= shrrr)- < 35 > 

The formula for Bernoulli polynomials is obtained by repeating exactly 
the same steps of the previous section. We will repeat these steps for the 
sake of clarity. 

The Bernoulli polynomials are usually defined using the generating func- 
tion 

te" ^ B t (x) 



nk=0 

or in terms of the function (f> x (t) as a generating function 



fc=i v 7 



Hence, by using the well-known identity B k (l — x) = (—l) k B k (x), we 
finally obtain 

(-l) k B k (x) = B k (l -x) = -^(O) = -I k -^ x (-t)] t=0 . (38) 
Since, for < t < oo, 

= E (^7 + — ) ^ - e_ 'r ^ ( 39 ) 
f Vn + 1 n / 

n=l 

uniformly, we may apply the fractional derivative operator termwise to 
obtain 



oo 

r s [M-t)] t=0 = J2 r i^Y + " et ) n_1 ^] t =o- (40) 

n=l 
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For —s = k — l, 



t\n—l xt~\ 



e lt=0 



n=l 

- ti^^^-er-^U (41) 



But 



Jfe-1 



lk 

^[(1 - eY^e-^ = 0, for n > k + 1 and 
(1 _ e'f-V 4 = J2(-l) J ( n ~ ^e^*; 

therefore, the infinite sum in f|5B|) reduces to a finite sum 



n — 1 o— n \ J / 

(j + x) 



n=l j=0 

n-l 

k-1 



n=l j=0 

^ \n(n+l) n 2 ) \j J 



j 

k 



*)* E (~T~r7\ + ^) A^(A0, (42) 
^— ^ \n(n + f ) n" 1 / 



where 



A ri ,,(fc) = ^(-l) J ( n )j(j+^-l)^ 1 - 
3 -=i w 



(43) 



The last equation combined with f)38p gives the explicit formula for 1^(1- 

x): 
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Bk(l -*)= ("l) fe+1 E {~r~~Tr\ + X —T-) A nAk), k > 1. (44) 

\n(n + 1 vr J 

n=l v y 

6. Bernoulli Numbers and Polynomials of Complex Index s 

There are many generalizations of integer-indexed Bernoulli numbers and 
polynomials to complex-indexed quantities. The reader may consult for ex- 
ample [l| or H and the references therein. Here, we approach the general- 
ization using fractional derivatives. 

The Bernoulli numbers and polynomials B s (x) for s complex can be de- 
fined using the fractional derivative operator of order 1 — s (i.e. replace k — 1 
by s — 1). When Re(s) > equation (1T5]) applies. Formulas dHJ and (J25J) 
yield the following natural definition when Re(s) > 0: 



t=o 



n + 1 

n=l 



oc 



n + 1 

n=l 

= (45) 
As for Bernoulli polynomials, their extension is obtained as follows 



t\n—Xxt] 

e Jt=0 



b.(i-*) = -/- 1 [^H)U = -E(^r+^r) /s_1 [(i-e 

n=l 

n=l j=0 V 7 



Jt=0 



oc 



E 5' n (l - s,x) 
n + 1 

n=l 

= <(l-s,x), (46) 

Thus, from equation (|46p . we see that the Bernoulli polynomials extend 
naturally to the entire function s((l — s,x). This is an illustration that entire 
functions are natural generalization of polynomials. 
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7. The Second Explicit Formula for B k 

In this section we will prove formula fl6]) using the globally convergent 
series representation of the Riemann zeta function discovered in 0. The 
series is given by 

n=0 

S n (s) being defined in f )23|) . 

Ii is easy to show that the sum S n +i(s) can be rewritten as 

S n+ i(s) = f^J (1 - e-ye-HT 1 dt, (48) 
and that for Re(s) > 0, 

1 f°° 

+ = — / ^f-'dt, (49) 
r 0) Jo 

where the function ip{t) is given by 

n=0 

Using the generating function ^(i), the Bernoulli numbers are now given 

by 

B k = V (k) (0) = -I k [v(-t)] t=o - (51) 

Again, we may apply the fractional derivative operator (— s = k) termwise 
to obtain 



n=0 

oo 



I k [(l - e') n e 
n + 1 



n=0 

CO 



y— 



jfc+i 



n=0 



n + l) 2 dt k+ 



iK 1 ~ e ) ] t =o 



(52) 
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But J3ft[(1 - e*)" +1 ] 4=0 = 0ifn+l>£; + 2. Therefore, the infinite sum 
in ( 1321 reduces to a finite sum 



n=0 

n=0 v ' j=0 

fc 1 n+1 / i 1 \ 

- (-^e^b-^C 1 ; 1 y + ' 

n=0 v ; j=0 \ J / 

k 

= (- l ) k+1 T,7^TV^ + i(k + l). (53) 



n " 



With an appropriate change of variable in the summation index, the last 
equation combined with ( 151|) gives the explicit formula for 5^: 



fc+i -, 

B k = (-l) fc+1 - 2 ^n{k + l),k> 0. (54) 

n=l 

The extension of the last explicit formula to Bernoulli polynomials and 
Bernoulli numbers of fractional index is straightforward. 



References 

[1] P. L. Butzer et al, Bernoulli Numberes and Polynomials of Arbitrary 
Complex indices, Applied mathematics letters, Vol. 5, No. 6, (1992), 
pp. 83-88. 

[2] K. Dilcher, Bernoulli Numbers Bibliography (1713- 
1990), Queen's Papers in Pure and Applied Mathemat- 
ics, No. 87, (1990). An online version is also available at 
http : / /www .mscs . dal . ca/~dilcher/berno ulli .html| (2007). 

[3] L. Fekih-Ahmed, On the zeros of the Riemann Zeta function, 
|http : //hal . archives-ouvertes . f r/hal-00476252_v3/[ (2011). 



12 



[4] L. Fekih- Ahmed, On the Hurwitz Zeta Function, 
http : //hal . archives-ouvertes . f r/hal-00602532[ (2011). 

[5] H. W. Gould, Explicit Formulas for Bernoulli Numbers, The American 
Mathematical Monthly, Vol. 79, No. 1, (1972), pp. 44-51. 

[6] H. Hasse, Ein Summierungsverfahren fur die Riemannsche Reihe, Math- 
ematische Zeitschrift, Vol. 32, pp. 456-465, (1930). 

[7] H. Laurent, Sur le calcul des derivees a indices quelconques, Nouvelles 
annales de mathematiques, serie 3, tome 3, (1884), pp. 240-252. 

[8] L. Saalschutz, Vorlesungen iiber die Bernoullischen zahlen, ihren zusam- 
menhang mit den secanten-coefficienten und ihre wichtigeren anwendun- 
gen, Springer, (1893). 

[9] N. M. Temme, Bernoulli polynomials old and new: Generalizations and 
asymptotics, CWI Quarterly, Vol. 8, (1995), pp. 47-66. 

[10] J. Worpitzky, Studien uber die Bernoullischen und Eulerschen Zahlen, 
J. Reine Angew. Math., vol. 94, (1883), pp. 203-232. 



13 



